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QCD Calculations of Radiative B Decays 

M. Misiak 

Institute of Theoretical Physics, Warsaw University, Warsaw, Poland 



The current status of B — » Xs-y decay rate calculations is summarized. Missing ingredients at the NNLO level are 
listed. The global normalization factor and non-perturbative effects are discussed. Arguments are presented that 
results for the cutoff-enhanced perturbative corrections have been misused in Ref. Il5l by applying them in the re- 
gion E-y E [1.0, 1.6] GeV, which means that the corresponding numerical effect on B (j3 — » X s *y, E 1 > 1.6 GeV) 
is unreliable. 
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1. Introduction 

The motivation for precision studies of radiative B 
decays is well known. First, they are sensitive to new 
physics loop effects that often arise at the same order 
in the electroweak couplings as the leading Standard 
Model (SM) contributions. Moreover, the inclusive 
B — > X s ^i rate is well approximated by the perturba- 
tively calculable radiative decay rate of the b-quark. 
The CLEO 0, BELLE and BABAR @ measure- 
ments have been combined by HFAG |4J to get 

B(B -> X s7 ) cxp = (3.52 ± 0.23 ± 0.09) x 1CT 4 (1) 

for E 1 > 1.6 GeV. The corresponding SM prediction 
that was published two years ago 1 reads [5| 

B(B ->AVt) sm = (3.15 ± 0.23) x 1CT 4 . (2) 

Its consistency with Eq. ([1]) provides strong constrains 
on many extensions of the SM (see, e.g., Ref. Q). 

Resummation of large logarithms (ot s lnM^/m^)" 
in the calculation of the decay rate is most conve- 
niently performed after decoupling the electroweak 
bosons and the top quark. In the resulting effective 
theory, the relevant flavour-changing weak interac- 
tions are given by a linear combination of dimension- 
five and -six operators 2 



Oi, 2 = (sr iC )(ci», 

03,4,5,6 = 

7 = 

8 = 



em b 

16tt 2 
gm b 

16tt 2 



^T a b R G 



(current-current 
operators) 

(four-quark 
penguin operators) 

(photonic dipole 
operator) 

(gluonic dipole /q\ 
operator) >■ ' 



One begins with perturbatively calculating their 
Wilson coefficients C, at the renormalization scale 
fio ~ (l%,m t ). Next, the Renormalization Group 
Equations (RGE) are used for the evolution of 
down to the scale jib ~ m,b/2. Finally, the operator 
on-shell matrix elements are calculated at /j, b . 



1 More recent contributions are discussed in Sees. [ 

2 The specific matrices T; and r". can be found in Ref. [Tj]. 



The Wilson coefficient RGE are governed by 
Anomalous Dimension Matrices (ADM's) that are de- 
rived from ultraviolet divergences in the Feynman di- 
agrams with operator insertions. Around 20000 four- 
loop diagrams with O1.2 insertions have been calcu- 
lated in Ref. [|[ to make the large logarithm resumma- 
tion complete up to O [a 2 (a s lnM^,/mj) 1, i.e. at 
the Next-to-Next-to-Leading-Order (NNLO) in QCD. 
Including such corrections is necessary to suppress the 
theoretical uncertainty in Eq. J2]) down to the level 
of the experimental error in Eq. {!]). The numerical 
effect of the four-loop ADM's on the branching ratio 
amounts to around —4% for fib = 2.5 GeV. At present, 
all the relevant Wilson coefficients Cj (//(,) are known 
at the NNLO [1,0. However, evaluation of the matrix 
elements at this order is still in progress — see Sec. 2] 



2. The global normalization factor 

In order to reduce parametric uncertainties stem- 
ming from the CKM angles as well as from the c- and 
6-quark masses, one writes the branching ratio as fol- 
lows El 



BIB 



Xr.ev) 



v t * s v tb 



A ^ [Pi E ) + N(E )} 



ttC 



(4) 



where a e 



a 



on shell 



and N(E ) denotes the 



non-perturbative correction (see Sec. [5]). The m c 
dependence of B — > X c ev is accounted for by 



C = 



Vub 



cb 



T(B -» X c ev) 
T(B^X u eD) 



(5) 



while P(Eq) is defined by the perturbative ratio 



r(& -> x sl ) 



\V cb /V ub \T(b ^ X u ei?) 



v t *v tb 



6a c 



P(E ). 



(6) 

The NNLO expression for the phase-space factor 
C ((5|) is a known function of m c /m b as well as non- 
perturbative parameters that affect N(Eq), too. All 
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these quantities are determined in a single fit from the 
measured spectrum of the inclusive semileptonic de- 
cay B — > X c ev. The fits are performed using either 
the IS or the kinetic renormalization schemes. The 
corresponding results for C and m c read 



C 

m c (m c ) 



0.582 ±0.016 
1.224 ±0.057 
<Ufi+ - 023 

U.O4O_ 033 



IS [11], 
kin. [12, 13]. 



1.267 

_ (7) 
The above MS-scheme values of m c have been ob- 
tained from the IS- and kinetic-scheme ones using the 
three-loop and two-loop relations, respectively. The 
three-loop relation for the kinetic scheme is not yet 
known. The ratio C is scheme-independent, but it is 
affected by the so-called weak annihilation contribu- 
tion By/A that remains unknown. Since By/A cancels 
out in Eq. Q , fixing its value is a matter of convention 
in the present context. Here, we follow the convention 
of Ref. [H}, namely _B\va(m = m b/2) = 0. 

The difference between the two determinations of 
C amounts to 1. 6a when counted in terms of the 
upper error of the very recent kinetic-scheme result 
[l2| . It is a consequence of using different experi- 
mental data sets, methodology and renormalization 
schemes. Fortunately, the effects of changing C and 
m c partially compensate each other in Eq. (|U) because 
d/dm c P{E ) < 0. For E = 1.6 GeV, I find 



B(B -> X sl ) = 



(3.15 ±0.23) x 10~ 4 , 
(3.25 ±0.24) x 10~ 4 , 



IS , 
kin. 



n (8) 

where the first result is just that of Ref. Q, while the 
second one has been obtained using the same code but 
with the input parameters from Rcfs. [la Il3j. 

The actual value of B (B -> X sl ) in Ref. [TJ is 
somewhat larger than 3.25 x 10~ 4 because P(Eq) was 
calculated there using the one-loop rather than two- 
loop determination of m c (m c ) from mjf n . In princi- 
ple, using the one-loop relation is allowed at 0{o? s ) 
because P{Eq) becomes m c -dependent only at 0{a 8 ). 



3. Cutoff-enhanced corrections 

The perturbative ratio P(E ) in Eq. © depends on 
the cutoff energy Eq via the dimensionless parameter 



2Eo 
mb 



(9) 



For very small <5, i.e. close to the kinematical endpoint 
Eq = mb/2, the usual ("fixed-order") perturbative ex- 
pansion breaks down because the corrections behave 
like powers of ln<5. In that region, one needs to re- 
sum large logarithms of S. Such a resummation of 
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Figure 1: Approximate cancellation of the logarithmic 
(4> < l' 1 ) and non- logarithmic ) terms in <f>^ away from 
the endpoint. The exact expressions from Eq. are 
represented by the solid lines. The Taylor expansions of 
and 4>^ around Eo = up to O(Eq) are shown by 
the dashed lines. 



the cutoff-enhanced corrections (i.e. corrections en- 
hanced by powers of In 5) has been performed up to 
the NNLO in Refs. [3, EE]. These results constitute 
a valuable contribution to our knowledge of the pho- 
ton energy spectrum in the endpoint region. However, 
they need to be treated with extreme care further from 
the endpoint, where logarithms of S no longer domi- 
nate. Naively, one might expect that resummation of 
small logarithms does not hurt, even if it is not an im- 
provement. Unfortunately, this is not the case because 
the logarithmic and non-logarithmic terms undergo a 
strong cancellation away from the endpoint. 

In order to illustrate this issue in a simple manner, 
let us consider only the dominant photonic dipole op- 
erator Oj in Eq. ©. When all the other operators 
are neglected, the fixed-order expression for the cutoff- 
dependence of P(Eq) is given by 



P(Eq) 
P(0) 



= 1 



(10) 



Each of the functions (ft^ can be split into two parts: 
(j)^ that is polynomial in In 5, and </>^ that vanishes 
at the endpoint. The explicit expressions for k = 1 
read [H 



} N ) 



^(S) = --\n 2 S--lnS--, 
3 3 9 ' 



W(6) = 



y*±^ 2 -^ 3 + ^-4)ln5. (11) 



In Fig. ED 



(f)j^ and their sum are plotted as 
functions of Eq. The endpoint is located at m&/2 ~ 

2.34 GeV. One can see that <fi^ begins to dominate 
around Eq — 2 GeV. On the other hand, already in 
the vicinity of Eq = 1.6 GeV, the cancellation of the 
two components of 4>^ is very strong. 
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Eq[ GeV] 

2.5 3 



Figure 2: Same as Fig. [T] but for 4> in the case when 
a s = a n J^ i (m h ) in Eq. (fTOJ. 



One may wonder whether similar cancellations oc- 
cur at higher orders, too. A positive answer concern- 
ing </>( 2 ) is immediate because this function can easily 
be derived from the results of Ref. fl7l |. The corre- 
sponding plot is presented in Fig. [2] for the case when 



a s = a s '~ (m b ) in Eq. |T 

No explicit results at order 0(a 3 ) are available. 
However, we know on general grounds that all the (f>^ 
behave like (2Po/mf,) 4 at small Eq. Two powers of Eq 
originate from P Mt/ in the vertex O-j in Eq. ([3]) , and two 
additional powers come from the phase-space measure 

(k) 

EydEy. Consequently, the Taylor expansions of <j> L 

and <j>N at small Eq up to 0{Eq ) must exactly cancel 
each other. These Taylor expansions are shown by the 
dashed lines in Figs. Q] and El One can see that both 
4>^ and (j)^ are well approximated by the dashed 
lines in the region below 1.6 GeV. It must also be the 
case at higher orders because (fr^ are polynomial in 
In 5 that is well approximated by the same expansion 
in the considered region (see Fig. [3]). Thus, cancella- 
tions like those shown in Figs. [T] and [2] are expected 
to occur at any order in a s (see "Note Added"). 
In the approach of Refs. [l4l fill ], logarithms of S 

(k) (k) 

have been resummed at the NNLO in (fr L , while 4> N 
have been retained in the fixed order. More precisely, 
Eq. (TTUj) has been effectively re-expressed as 



P(Eo) 
P(0) 



X 



6 (1) 



, (2) 
'N ' 



(12) 



and X has been calculated up to O (a 2 , x a™ ln m 6) , 
with n — 0,1,2,..., and m = n, n + 1, . . . , m max (n). 
This is a reasonable approximation only in the re- 
gion very close to the endpoint where no cancella- 
tion between the two components takes place. Else- 
where, it leads to overestimating the numerical ef- 
fect of the 0(a^) terms in Eq. (llOj) by a factor of 



order <^ 3) /0 (3) 

around 4.6 for 
Eq = 1.0 GeV. 



En 



• (mb/(2Eo)) that amounts to 
= 1.6 GeV, and around 30 for 
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Figure 3: In S (solid) and its expansion (dashed) up to 
O(E 3 ). 



Unfortunately, it was precisely the range 
E 1 G [1.0, 1.6] GeV where the authors of Ref. [l|| 
applied their results to calculate the effect on 
B{B -»X S 7, £ 7 > 1.6 GeV). They adopted the 
fixed-order result at Eq = 1.0 GeV from Ref. [j| 



B(B X sl , E y > 1.0 GeV) = 3.27 x 10" 



(13) 



and supplemented it with their own numerical value 
of [P(1.0) - P(1.6)] /P(Q) that follows from Eq. (p)). 
That value is almost twice larger than in the fixed- 
order NNLO calculation. In the end, their result for 
the branching ratio with a cutoff at Eq — 1.6 GeV 
was considerably lower than the one in Eq. ([2]). In 
view of the above remarks, their prediction should be 
considered unreliable. 



4. Status of the NNLO QCD calculations 
of the matrix elements 

On the l.h.s. of Eq. © that defines P{EnV the 
denominator is already known at the NNLO [ill [l^| . 
In the expression for the numerator 



T{b -> X sl ) E „ i>Eo 



G 2 F a em ml 
32tH 



\v tb v; s { 



x £ Cf Cf G^EqUIA) 

the quantities Gy are determined by the matrix ele- 
ments of Ox, . . . , 0$. 

So far, only G77 has been evaluated up to 0{a 2 s ) in 
a complete manner [13, HO, The remaining 
are fully known at the Next-to-Leading Order (NLO), 
i.e. up to 0(a s ) (see Ref. 22] for a complete list of 
references). At the NNLO, it is practically sufficient 
to restrict considerations to Gij with i,j S {1,2,7, 8} 
because the four-quark penguin operators have small 
Wilson coefficients. It is often convenient to apply 
the optical theorem, and calculate Gy by summing 
imaginary parts of the &-quark propagator Feynman 
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Table I Additional known effects not included in Eq. JB). 



Figure 4: Examples of Feynman diagrams that contribute 
to G77, G27 and G22 at the NNLO. The dashed vertical 
lines mark the unitarity cuts. Black squares represent the 
operators O2 and O7 in Eq. ©. 



diagrams. One can see in Fig. [5] that imaginary parts 
of three-, four- and five- loop diagrams occur at the 
NNLO in G77, G27 and G22, respectively. 

A relatively simple set of the NNLO contribu- 
tions is given by diagrams with quark loops on the 
gluon lines. The quark in the loop is either mas- 
sive (charm and bottom) or treated as massless (up, 
down and stran ge) . Such contributions are already 
known 0, HI HI, HI for all the with i,j £ 
{1,2,7,8}, except for the massless case in Gis and 
(?28- The BLM [26| (or large- (3q) approximation for 
the complete NNLO correction is derived from the 
massless quark results [H, [24[ ■ 

In Ref. [27|, the asymptotic behaviour for m c ^> 
mt,/2 was calculated for all the non-BLM NNLO cor- 
rections to Gij with i,j £ {1,2,7,8}, except for G78 
and Ggs- Next, an interpolation in m c of these correc- 
tions was performed assuming that the interpolated 
quantities vanish at m c — 0. The result of that pro- 
cedure was an essential input for the NNLO estimate 
in Eq. ([2]). The overall error of around 7% in the 
branching ratio was obtained by combining in quadra- 
ture four types of uncertainties: 5% non-perturbative, 
3% parametric, 3% higher-order, and 3% due to the 
m c -interpolation ambiguity. 

The results in Eqs. and © do not include several 
contributions to the branching ratio that are known 
at present. These additional effects are summarized in 
Tab. [I] They sum up to around +1.6%, which is small 
when compared to the overall uncertainty of around 
7%. Therefore, Eq. (jHJ) can still be treated as an up- 
to-date SM prediction. 

As the reader has already noticed, even in the mc- 
interpolation approach of Ref. [27| , there are still some 
missing NNLO ingredients, namely: 

• the BLM contributions to Gis and G28, 

• the large-m c results for G78 and Gss- 

Their numerical effect on the branching ratio is ex- 
pected to remain within the estimated higher-order 



The BLM terms from Ref. [24] 

O s in the 4-loop ADM's [8] 

b and c loops on gluon lines [19 . 21 . 251 

Non-perturbative Q(a s A/mb) effects [26 

Non-perturbative collinear effects [291 



Total 



+2.0% 
-0.3% 
+1.6% 
-1.5% 
-0.2% 



+1.6% 



uncertainty of around 3%. The calculation of the 
most interesting G78 is very advanced [3(ifl . and the 
results should become available soon (for any value of 
m c ). As far as Gss is concerned, its calculation in the 
large-m c limit will automatically give the result for 
any value of m c . 

For the full NNLO calculation, the currently miss- 
ing ingredients (apart from the ones listed above) are 
the non-BLM corrections to 

(i) G17 and G27, 

(ii) Gn, G12 and G22, 

(Hi) Gis and G28- 

The calculation of (i) in the m c = limit is quite 
advanced, but also extremely difficult and time- 
consuming — see Ref. [3l| for the status reports. Once 
it is finished, the main challenge will amount to find- 
ing (ii), even for m c = 0. The corrections (Hi) are 
expected to be numerically less important. 

When the non-BLM corrections are known at 
m c — sometime in the future, the interpolation in 
m c is still going to be necessary. However, our error 
estimates should become more solid once the BLM 
approximation is no longer used at the boundary. 

Finding all the non-BLM corrections for the actual 
value of m c ~ m^/A is even more difficult, but it must 
be considered at some point, too. A calculation of 
the IR-divergent two-particle-cut contributions to (i) 
is being currently performed for arbitrary m c [311 ]. 
The IR-convergent two-particle-cut contributions to 
(ii) for arbitrary m c are already known because they 
are given by products of the NLO corrections. The 
(n > 3)-particle-cut contributions to (ii) vanish at 
the endpoint, so their numerical relevance should be 
diminished by the high photon energy cutoff. 

Apart from the NNLO corrections, there are other 
perturbative contributions to T(b — > X s j)e~>e that 
have been neglected so far, namely tree-level dia- 
grams with the u-quark analogues of O1.2 and the 
four-quark penguin operators 03,4,5.6- Such contri- 
butions are suppressed with respect to the leading 
term by |(C? >2 , C 3A ^ e )/C 7 \ 2 < 0.2, where G" 2 = 
{V* s V ub )/(V t *Vtb)C 1: 2, as well as by the high photon 
energy cutoff. A quantitative verification of how small 
they really are should become available soon 32 1. 
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5. Non-perturbative effects 

Let us begin with considering a simplified world 
where m c — mj. There, in the decay of the B meson, 
a high-energy photon (I? 7 ~ mb/2) can be produced 
in four different ways: 

1. Hard: The photon is emitted directly from the 
hard process of the &-quark decay. 

2. Conversion: The 6-quark decays in a hard way 
into quarks and gluons only. Next, one of the 
decay products scatters in a non-soft radiative 
manner with the remnants of the B meson. This 
can be viewed as a parton-to-photon conversion 
in the QCD medium. 

3. Collinear: In the process of hadronization, a 
collinear photon is emitted. 

4. Annihilation: An energetic qq state produced in 
the B meson decay disintegrates radiatively. 

The hard way would be the only one if no other 
operators but O7 were present in the effective theory. 
Non-perturbative effects in such a case were first an- 
alyzed in Ref . [33] . They arise as corrections of order 
A 2 /m\ to T(b -> X s ~f) El> E when (m b - 2E ) ~ m b . 
Moreover, these 0(A 2 /m 2 ) terms cancel out in N(Eq) 
in Eq. Q with the analogous non-perturbative cor- 
rections to the charmless semileptonic rate. Thus, we 
are left with the small C(A 3 /mg) effects The 
O (A 2 / (m& — 2£'o) 2 ) corrections are also small [35[ for 
Eq < 1.6 GeV. All these terms are included in Eq. ([8]), 
and affect the branching ratio by around —0.7%. 

The photon production via conversion is suppressed 
both by a s (due to the non-soft scattering) and by 
A/ nib (due to dilution of the target). The analysis 
in Ref. [28| shows that no other suppression factors 
occur. An effect on the branching ratio of roughly 
— 1.5 ± 1.5% was found in that paper (see Tab. Q}. 

In our simplified case (m c = to;,), the tree-level de- 
cay b — > (quarks & gluons) s is possible only via the 
operators 03,4,5.6,81 an d by the u-quark analogues of 
Oi,2- Consequently, the collinear photon emission 
is suppressed either by a s |C 8 /CV| 2 < 0.08, or by 
|(C" 2 , C 3 ,4, 5 , 6 )/C7| 2 < 0.2. Moreover, there is an ad- 
ditional suppression by products of the quark electric 
charges and, most importantly, by the high photon 
energy cutoff. The results of Ref. [29[ lead to an es- 
timate that the non-perturbative collinear effects due 
to C*8 amount to around —0.2% in the branching ratio 
for E = 1.6 GeV (see Tab. 0J. 

As far as annihilation is concerned, photons orig- 
inating from decays of 7r°, 77, rj' and ui are removed 
on the experimental side as (huge) background. Con- 
tributions from the other established qq mesons are 
negligible. The corresponding perturbative diagrams 
are responsible for only around 0.1% of the total rate 
for E = 1.6 GeV. 



Once the assumption m c — nib is relaxed, the 
C(A 2 /m 2 ) hard effects from Oi,2 get replaced by a 
series of the form [36] 




with quickly decreasing coefficients b n . The calculable 
leading term has been included in Eq. (|8|). It affects 
the branching ratio by around +3.1%. 

All the quantitatively estimated non-perturbative 
effects that have been mentioned so far sum up to 

(-0.7 - 0.2 - 1.5 + 3.1)% = +0.7%. (16) 

However, since their evaluation is often very uncer- 
tain, and the knowledge of 0(a s A/rrib) contributions 
is by no means complete, a non-perturbative error of 
±5% has been assumed in Eq. ([8]), as already men- 
tioned Sec.Q] Probably the most interesting of all the 
unknown 0(a s A/mb) effects originate from charm an- 
nihilation in the massive (cs)(qc) intermediate states 
(g = u or d), for the actual value of m c ~ nib / 4. 

One should remember that the error in Eq. ([T]) 
is affected by a non-perturbative theoretical uncer- 
tainty, too. It follows from the fact that the actual 
measurements are not performed with Eq ~ 1.6 GeV. 
The most precise experimental results correspond to 
significantly higher photon energy cutoffs for which 
the 0(A n /(nib — 2E ) n ) effects are no longer small. 
These effects are described by a non-perturbative 
shape function [53] that is constrained by the semilep- 
tonic data. The very recent analysis [381 ] of this func- 
tion and its effects on the B — » X s -f photon spectrum 
can hopefully provide input for the future HFAG av- 
erages. 



6. Conclusions 

Thanks to the efforts of the past years, the un- 
certainties in B(B — > A s 7j have reached the level of 
around ±7% on both the experimental and theoretical 
sides. A significant progress in the perturbative cal- 
culations is expected in the near future. However, un- 
derstanding the 0{a s A/nib) non-perturbative effects 
remains the key issue. 



7. Note Added 

After the first version of this article was submitted 
to the arXiv, Einan Gardi pointed out to me that 
the approximate cancellation of logarithmic and non- 
logarithmic terms (Sec. [3]) has already been discussed 
in Ref. 
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